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ABSTRACT 

The GGD27 complex includes the HH 80-81-80N system, which is one of the 
most powerful molecular outflows associated with a high mass star-forming region 
observed up to now. This outflow is powered by the star associated with the 
source IRAS 18162-2048. Here we report the detection of continuum emission 
at sub-arcsec/arcsec resolution with the Submillimeter Array at 1.36 mm and 
456 /im, respectively. We detected dust emission arising from two compact cores, 
MM1 and MM2, separated by about 7" (-12000 AU in projected distance). MM1 
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spatially coincides with the powerful thermal radio continuum jet that powers 
the very extended molecular outflow, while MM2 is associated with the protostar 
that drives the compact molecular outflow recently found in this region. 

High angular resolution obervations at 1.36 mm show that MM1 is unresolved 
and that MM2 splits into two subcomponents separated by ~ 1". The mass of 
MM1 is about 4 M and it has a size of < 300 AU. This is consistent with MM1 
being associated with a massive and dense (n(H 2 ) > 10 9 cm -3 ) circumstellar 
dusty disk surrounding a high-mass protostar, which has not developed yet a 
compact HII region. On the other hand, the masses of the two separate com- 
ponents of MM2 are about 2 M each. One of these components is a compact 
core with an intermediate-mass young protostar inside and the other component 
is probably a pre-stellar core. 

MM1 is the brigthest source at 1.36 mm, while MM2 dominates the emission 
at 456 jam. These are the only (sub) millimeter sources detected in the SMA 
observations. Hence, it seems that both sources may contribute significantly to 
the bolometric luminosity of the region. Finally, we argue that the characteristics 
of these two sources indicate that MM2 is probably in an earlier evolutionary 
stage than MM1. 

Subject headings: circumstellar matter — ISM: individual (GGD27, HH 80-81, 
IRAS 18162-2048) — stars: formation — stars: early type — submillimeter: ISM 

1. INTRODUCTION 

It is well known that low-mass stars form by accretion via circumstellar disks. In 
addition, collimated outflows are believed to remove excess angular momentum from the 
system. However, the formation process of high-mass stars (> 8M ) is yet unclear (see 
Zinnecker & Yorke 2007, for a recent review). When a protostar accretes enough mass to 
ignite hydrogen, it starts emitting a large amount of ionizing photons. This radiation ionizes 
the surrounding gas thereby forming an H II region. At first glance, the pressure from 
the stellar radiation and the ionized gas, could prevent further accretion (see e.g., Larson 
& Starrfield 1971; Kahn 1974). However, accretion through disks may circumvent these 
problems, allowing protostars to grow up to 60-100 M (Nakano 1989; Jijina & Adams 1996; 
Yorke & Bodenheimer 1999). Several theoretical models have been proposed to explain 
the formation of high-mass stars (see Zinnecker & Yorke 2007 for a review). Models based 
on accretion require high accretion rates (e.g., McKee & Tan 2003), and/or accretion of 
ionized material after the hydrogen burning has started in the star (Keto & Wood 2006). 
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More dynamical models are based on the interaction between the collapsing fragments of 
the parental cloud (i.e., competitive accretion, Bonnell et al. 2001), or coalescence of less 
massive protostars (Bonnell, Bate, & Zinnecker 1998). At present, the observational data 
do not allow to distinguish between these models. 

Because of their short lifetimes, high-mass protostars remain embedded in their natal 
envelopes. This situation hampers the detection and characterization of disks, since disk 
emission is usually mixed with that of its own envelope. In addition, high-mass protostars 
are located typically a few kiloparsec away and they usually form in clusters, which makes 
it difficult to resolve without sufficient angular resolution. Hence, currently there are only a 
handful of appropriate candidates for detailed studies. 

There is some observational evidence of disks around B-type protostars (with stellar 
masses < 20 M , see e.g., Cesaroni et al. 2007; Zhang 2005). Disks around O-type protostars, 
are elusive. Instead of disks, O-type protostars (with masses > 20 M ) are surrounded 
by ~ 10000 AU rotating molecular structures, which seem to be dynamically unstable or 
collapsing structures (e.g., Beltran et al. 2006; Zapata et al. 2009). 

In the Sagittarius region, at a distance of 1.7 kpc (Rodriguez et al. 1980), one of the 
most powerful radio continuum jets has been associated with high mass star formation. The 
HH 80-81-80N system, as part of the GGD27 complex, has several remarkable characteristics: 
(i) the length of the thermal jet, ~5.3 pc, is one of the largest found so far (Marti, Rodriguez, 
& Reipurth 1993). Multiwavelength VLA studies (Marti et al. 1993; Marti, Rodriguez, & 
Reipurth 1995, 1998, 1999) show that this radio jet ends at two Herbig-Haro objects to the 
south (HH80 and 81; Reipurth & Graham 1988) and at a radio source to the north (HH 
80 North; Marti et al. 1993; Girart et al. 1994, 2001; Masque et al. 2009). (ii) The HH 
80-81 objects are very bright in the optical (Heathcote, Reipurth, & Raga 1998) and radio 
wavelengths (Rodriguez & Reipurth 1989), ~ 10 times brighter than those of classical HH 
objects (Curiel 1995). (iii) The motions found in these HH objects and related ejections are 
very large (optical line widths up to 700 km s -1 , Reipurth & Graham 1988; proper motions 
up to ~500 km s _1 , Marti et al. 1998). 

This radio jet may be driving a molecular outflow detected by single dish CO observa- 
tions (Yamashita et al. 1989; Ridge & Moore 2001; Benedettini et al. 2004). The northern 
lobe of the outflow approaches us (blue lobe) while the southern one moves away from us 
(red lobe). Recent Spitzer observations show a V-shaped biconical structure at 8 /xm roughly 
coincident with the center of the molecular outflow (Qiu et al. 2008). The vertex of this 
biconical structure coincides with the radio jet. 

The powering source of the radio jet and the CO outflow, coincides with the IRAS source 
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18162-2048, which has been observed at 20 cm (Marti et al. 1993) and at 5.8 /im (Qiu et al. 
2008). Based on its IRAS fluxes, the bolometric luminosity of this source is ~ 1.9 x 10 4 L , 
corresponding to a BO zero age main sequence (ZAMS) star. Other typical tracers of young 
massive stars have been found in this region, such as compact CS emission (Yamashita et al. 
1991) and bright H 2 and CH 3 OH masers (Rodriguez et al. 1978, 1980; Gomez, Rodriguez, 
& Marti 1995; Val'tts et al. 2000; Kurtz, Hofner, & Alvarez 2004). Ground-based infrared 
observations have identified this source as part of a cluster of intermediate-to-high mass 
young stars (Aspin & Geballe 1992; Aspin et al. 1994; Stecklum et al. 1997) and at the 
same time as the responsible for illuminating the reflection nebula seen at these wavelengths 
(Aspin et al. 1991; Aspin & Geballe 1992). 

BIMA observations at 1.4 and 3.5 mm show compact dust emission (< 3" in size), as 
well as compact SO(5s — 4 4 ) line emission (< 5" in size) toward the driving source of the 
outflow (Gomez et al. 2003). These observations also show that there is a second millimeter 
source about 7" to the NE from the center of the radio jet, with weak continuum emission 
at 1.4 mm, coinciding with warm NH3(1,1) and HCN(l-O) emission. This source nearly 
coincides with a weak radio continuum compact source detected at 3.5 cm and a strong 
water maser spot (named VLA3 in Gomez et al. 1995). It is also associated with a very 
young one-side molecular outflow observed at high velocity towards the south-east (Qiu & 
Zhang 2009). 

In this paper, we report Submillimeter Arra}0 (SMA; Ho, Moran, & Lo 2004) observa- 
tions of the GGD27 system at 1.36 mm and 456 /im. Pushing at longer wavelengths allows 
the study of colder (10-100 K) regions. The SMA has the ability of observing at two fre- 
quencies simultaneously, sampling the (sub)millimeter spectrum of the sources with a single 
experiment. The SMA also allows higher angular resolution (< 1"), which is crucial for defin- 
ing structures and resolving the multi-component systems. The SMA also has a very wide 
spectral band and can detect many spectral lines simultaneously. We have applied different 
phase calibration techniques to the submillimeter data. We also obtained spectral-line data 
which we will publish elsewhere. The observations are described in §2 and the results are 
presented in §3. A discussion on the analyzed data is presented in §4 and finally, a summary 
with the main conclusions is presented in §5. 



lr The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. OBSERVATIONS 

The observations were conducted in two epochs separated by two years (2005-2007). 
The details of the observations are shown in Table [fl During the first epoch we ob- 
served with the SMA in its compact configuration while using the dual frequency mode 
(220/658 GHz). During the second epoch we observed only at 220 GHz with the SMA in 
its very extended configuration. The phase center was RA(J2000.0)= 18 fc 19 m 12!099 and 
DEC(J2000.0)= — 20°47'30 / .'00, near the position of the source at the center of the thermal 
radio continuum jet of the GGD27 system. 

On 2005 August 24, the two SMA receiver bands were tuned simultaneously at 215/225 GHz 
and 648/658 GHz (each receiver band has two sidebands with a bandwidth of 1.968 GHz 
each, and separated by 10 GHz; they are labeled lower and upper sideband respectively, 
LSB/USB). The interferometer was in its compact configuration and the projected baselines 
ranged from 13 to 68 m. This configuration produced a synthesized beam size of 8'.'1 x 3"0 
at 220 GHz and of 2'!4 x l'.'l at 658 GHz. The SMA primary beam has a size of about 50" 
and 17" at 220 GHz and 658 GHz, respectively. Weather was good during the observations, 
with a relative humidity of 7% and a zenith opacity ranging between T225 = 0.05 and 0.08 
(t"690 ~0.9 - 1.5, following Masson 1994), measured with the radiometer of the nearby Caltech 
Submillimeter Observatory, at 225 GHz. The recorded data were edited, calibrated, imaged 
and analyzed using the Miriad (Sault, Teuben, & Wright 1995) and AIPS (developed by 
NRAO) packages. 

On 2007 May 29, we observed with the low frequency receiver tuned at the same fre- 
quency as in the previous observation (215/225 GHz). The array was in its very extended 
configuration with projected baselines ranging between 50 and 514 m, yielding an angu- 
lar resolution of 0"7 x 0"4. At this time, the zenith opacity during the observations was 
T225 = 0.19 on average. This second set of data was also edited, calibrated, imaged and 
analyzed using the Miriad and AIPS packages. 

To image the data we applied a u-v weighting during the Fourier transformation stage. 
The robust parameter (Briggs 1995; Briggs, Schwab, & Sramek 1999) controls this weighting, 
relating the weight of a u-v cell to the number of visibilities within it. The robust parameter 
ranges between -2 and 2 in Miriad. A robust of 2 corresponds to a natural weighting and a 
robust of -2 corresponds to a uniform weighting. While the natural weighting gives a lower 
rms noise level, the uniform weighting minimizes the sidelobe levels and gives a better spatial 
resolution. Values of the Miriad's robust parameter between -2 and 2 provide a compromise 
between these two extrema. 



-6 - 



2.1. SMA 220 GHz observations 

2.1.1. Compact configuration 

At 220 GHz the system temperature was between 70 and 150 K, mainly depending on 
the elevation of each source along the ~9 hours observing track (about 3 hours of which were 
on target source). The quasar 3C454.3 was used to correct the bandpass and Callisto was 
used to do the absolute flux calibration. The observations were carried out in cycles of 25 
minutes, observing VX Sgr for 10 minutes (in order to try a phase-transfer calibration), then 
moving to GGD27 for 10 minutes and finally observing the quasar J1911-201 for 5 minutes 
in each cycle. The gain solutions obtained for the quasar J191 1-201 (with flux densities 
of 2.13 and 2.18 Jy at 215 and 225 GHz, respectively) were applied to GGD27 and VX 
Sgr. The flux calibration was performed with the visibility amplitudes from Callisto and its 
uncertainty was estimated to be less than 20% by comparing measured and expected fluxes 
for several observed planets and moons (Ganymede, Neptune, Ceres and Callisto; see Table 
[2]). Dirty maps were obtained after discarding one of the seven initial antennas because it 
was unable to follow the atmospheric-phase behavior properly (probably due to instrumental 
problems). The channels with line emission were removed from the bandwidth before making 
a continuum channel for each sideband. Then, both sidebands (with similar rms noise levels) 
were combined and self calibration applied with an averaging time of 3 minutes. An image of 
the region was then cleaned and restored, applying natural weighting. The resultant image 
has a synthesized beam size of 8'.'1 x 3'.'0 (P.A.= 34? 1) and an rms of ~ 10 mJy beam -1 , 
which appears to be limited by the dynamic range. The astrometric accuracy was O'/l. In 
addition, nine molecular lines (five SO2 transitions, two lines of SO, one of C 17 and another 
of H2CO) were detected through both sidebands. The results of the analysis of these lines 
will be published in a subsequent paper. 

2.1.2. Very extended configuration 

During the 2007 May 29 observations, the system temperature ranged between 150 and 
400 K. 3C454.3 was used for bandpass calibration. The frequency setup was identical to 
the compact configuration and seven of the SMA antennas were used on this occasion. The 
phase and amplitude calibrator was also J1911-201, with a measured flux density of 1.24 Jy 
at 215 GHz and 1.16 Jy at 225 GHz. The fluxes at 225 GHz of 3C454.3 and MWC349 (which 
was also included as a flux calibrator) were 8.57 and 1.97 Jy respectively. Thus, comparing 
these values with those found for these calibrators in the SMA data base, the uncertainty 
in the flux scale was estimated to be less than 15%. After the visibilities were calibrated 
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and channels with line emission were removed, dirty maps of the continuum emission were 
constructed. Then the two sidebands were combined and self calibrated with an averaging 
time of 3 minutes. An image was cleaned and restored applying natural weighting. The 
synthesized beam of the resulting image has a size of CC/68 x (X'38 (P.A.= 13?0) and the 
rms of the image is ~3.0 mJy beam" 1 . The accuracy on the position of the sources based 
on phase noise was ~ Of! 02, which coincides with the absolute positional accuracies of the 
calibrators observed. At least seven of the nine molecular lines observed in the compact 
configuration were also detected at this epoch. 



2.2. SMA 658 GHz observation 

On 2005 August 24 we observed simultaneously with two receiver bands tuned at 
215/225 GHz and 648/658 GHz. At 658 GHz, the system temperature varied between 
2000 and 4000 K, more than one order of magnitude higher than that at 220 GHz due to the 
higher atmospheric emission at 658 GHz. The bandpass was calibrated finding a solution 
for each baseline of the quasar 3C454.3 (with the Miriad task UVGAINS), which had a flux 
density of 13.5 Jy at this frequency, while the amplitude gains were obtained with visibility 
data from Ceres, which had the best signal-to-noise ratio among all calibrators observed at 
this frequency. The estimated flux uncertainty is less than 15%. Astrometric accuracy is 
estimated to be 0".2. Given that interferometric observations at high frequencies, such as 
650 GHz, are still very difficult, observations were planned to be able to phase-calibrate 
the data: a) using the standard procedure (using the quasar J191 1-201 as the phase cali- 
brator), b) using the water maser at 658 GHz associated with a nearby evolved star (VX 
Sgr) to phase-calibrate the data, and c) applying a phase transfer from the simultaneous 
observations at 220 GHz. 



2.2.1. Standard calibration 

One of the main problems concerning submillimeter calibrations is the lack of strong 
calibrators at high frequencies. At the time of observations, the quasar J1911-201, used 
without problems as phase calibrator at 220 GHz, was too weak (< 1.5 Jy) to be detected 
at 658 GHz. Thus, the standard procedure of phase calibration at high frequency was not 
possible. 
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2.2.2. Calibration using a close by H2O maser 

The strong 658.007 GHz H 2 li,o-lo,i maser line of VX Sgr (a red supergiant) was used 
as phase calibrator for the USB high frequency sideband. VX Sgr is ~3° away from GGD27 
and this fact together with the fairly good weather, made the calibration cycle time of 25 
minutes good enough to remove the slow instrumental and atmospheric drifts. To find the 
best gain solutions, only the strongest channels with a stable phase were used (averaging 
the line channels into a continuum channel). The visibilities from 6 antennas were used 
to self calibrate the water maser (one of the antennas did not follow the phase calibration 
properly, probably due to instrumental problems). Then, the phase-solutions were applied 
to the GGD27 upper sideband data. A cleaned image of the upper sideband shows emission 
consistent with the sources detected at 220 GHz. Then, the continuum emission was self 
calibrated to correct for the rapid tropospheric fluctuations, improving the signal-to-noise 
ratio (SNR) by a factor of ~ 4. No spectral lines were detected at this frequency. The image 
was natural weighted, giving a synthesized beam of 2'.'41 x 1"14 with a P.A.= 37?1 (see Fig. 

Two compact continuum sources were clearly detected at this frequency. In order to 
properly determine the flux of the sources (see Table [3]), a correction for the attenuation of 
the primary beam was applied to the final image. The rms (measured into the primary beam 
and out of the sources) before and after the primary beam correction was 457 and 678 mJy 
beam -1 . 



2.2.3. Phase transfer from low frequencies 

A different way to do phase calibration was tried by transferring the phase solutions 
from 215 to 658 GHz. This method requires very good quality data and the stability of the 
response of the detectors. The phase calibration was tried using two strong masers of VX 
Sgr (the water maser at 658 GHz and a SiO maser at 215 GHz), but a high rms noise level 
of the phase at high frequency on GGD 27, prevented the success of the calibration of this 
source. However, we report the success of this calibration technique on the strong water 
maser of VX Sgr. The details of this calibration are presented in the appendix. 

2.2.4- Positional corrections 

Comparing the continuum images of GGD27, obtained at 220 and 658 GHz, we found 
that the peak positions of the two observed sources in the high frequency image were offset 
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about 0'.'7 to the NW with respect to their positions found in the low frequency image. In 
addition, the coordinates of the VX Sgr evolved star and its phase center did not match, 
showing a similar trend as the offset in the GGD27 image. Since we used the water maser 
at 658 GHz associated with VX Sgr for the phase calibration of GGD27, an artificial offset 
was transferred to the GGD27 image. 

In order to correct this positional offset, we obtained the peak position of the centroid 
of the 215 GHz SiO maser of VX Sgr and compared it with the coordinates of the phase 
center used to observe this evolved star. This comparison indicates that the SiO maser is 
located about 0'.'80 to the north and O'.'Ol to the west from the coordinates used during the 
observations. A similar shift is found when comparing the coordinates of this star from 
SIMBAD database and the coordinates of the phase center for VX Sgr (0'.'69 to the north 
and O'.'Ol to the west). Since the spatial distribution of the SiO and H 2 maser spots is small 
(< 0'.'4) compared with the synthesized beam size of the observations, we assume that the 
centroid of both types of masers coincide spatially. We applied the estimated positional offset 
of the centroid of the SiO maser to the calibrated water maser position at 658 GHz, obtaining 
that the new position of the H 2 maser coincides with the SiO maser position within the 
pixel size (0'.'18) at 658 GHz. Then, we applied the same correction to the GGD27 data. 
The coordinates of the two submillimeter sources in the corrected high frequency image 
spatially coincide with those obtained at 1.36 mm, within 0'.'2, the astrometric accuracy of 
the 658 GHz image. The high frequency image presented here (see Fig. [2]) was corrected as 
described above. 



2.3. Notes on flux calibration 

Table[2]shows the measured and expected flux densities for all the calibrators used during 
these observations. We find that in all epochs, the obtained flux densities coincide with those 
expected within a 15% error. The exceptions are Ceres for the 2005 August observations 
at 220 GHz, and Neptune for the 2005 August observations at 658 GHz (the latter was 
partially resolved by the interferometer). For the May 2007 observations towards 3C454.3 
(the strongest calibrator) we obtain an exceptional 3.8% accuracy at 220 GHz, compared to 
the interpolated data measured between 215 and 225 GHz from the SMA quasar monitoring. 
Finally, positional accuracy of the calibrators is better than O'.'l. 
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3. RESULTS AND ANALYSIS 

3.1. Continuum emission 

Figure [T] presents the 1.36 mm (220 GHz) continuum emission images towards GGD27 
obtained with the SMA in its compact and very extended configurations. The low angular 
resolution image shows a strong and elongated source, with a main peak and a possible second 
peak towards the NE. The high angular resolution image clearly separates the continuum 
emission into two sources, a strong and compact one, MM1, located near the phase center 
and a weak and extended one, MM2, located about 7" to the NE from MM1. The 456 
(658 GHz) image (see Fig. [2]) also shows two compact sources spatially coincident with MM1 
and MM2 (see Table [3]). MM1 dominates the emission at 1.36 mm and coincides spatially 
with the powerful radio continuum jet found in this region (Marti et al. 1993, 1995; Gomez 
et al. 1995). MM2 spatially coincides (within O'.'l) with a weak radio continuum source and 
a water maser (Gomez et al. 1995; Marti, Rodriguez, & Torrefies 1999; Kurtz & Hofner 
2005). 

Table E] shows the results of the Gaussian fits (position, peak flux, flux density and 
deconvolved size) carried out with the IMFIT task of AIPS in the three images, for MM1 
and MM2. The image at 1.36 mm obtained with the compact configuration (left panel of 
Fig. [1]), shows a single extended source, which was fitted simultaneously with two Gaussians. 
This fit gives fluxes of 610±8 mJy and 443±9 mJy for MM1 and MM2, respectively and 
deconvolved sizes of 3'.'5± O'.'l x 1'.'3±0'.'2; P.A.=176°±3° for MM1 and 4'.'9±0'.'2 x 1'.'5±0'.'5; 
P.A.=164° ± 3° for MM2. In the very extended configuration image (right panel of Fig. [T]), 
MM1 has a flux density of 441 ± 6 mJy and appears unresolved with an upper limit to its 
diameter of 0'.'36 (see Table |3] and section § 4 for the derivation of a more stringent upper 
limit). MM2 splits into two main components separated by ~ 1". The eastern component, 
MM2(E), is compact and has a flux density of 57±6 mJy, with an upper limit to its diameter 
of 0'.'36. The western component, MM2(W), is slightly weaker (50±10 mJy) and extended, 
with a deconvolved size of l'.'l ± 0'.'2 x 0'.'5 ± O'.'l; P.A.=34° ± 11°. In the 456 /im image (Fig. 
|2]) the detected sources are unresolved (diameter smaller than 1'.'7), with fluxes of 5.3±1.7 Jy 
for MM1 and 19±1.5 Jy for MM2. 

The total flux measured in the low angular resolution image at 1.36 mm (1.20 Jy) is 
slightly greater (by ~ 0.15 Jy) than the sum of the fluxes of MM1 and MM2 obtained 
by fitting two Gaussians. This could suggest the existence of a weak extended component 
surrounding both sources. In addition, the total flux is greater than that measured in 
the high angular resolution image (550 mJy), which indicates that the extended emission 
in the low angular resolution image is resolved out in the high angular resolution image. 
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In order to analyze the nature of this extended emission, the data at 1.36 mm of both 
configurations were combined (Fig. [3]). The low and high angular resolution data overlaps 
in a uv-plane region between 18 and 52 kA, in which their fluxes concide within 3%. This 
allowed a straightforward combination of the data, using the INVERT routine of Miriad; 
the resultant images were obtained using several robust parameters (2, 0.5, 0.25 and 0) to 
sample intermediate angular resolutions between the images of compact and very extended 
configurations. In the robust=2 image (see Fig. [3]), the emissions from MM1 and MM2 are 
blended due to the synthesized beam size. However, when angular resolution improves, the 
emission from each source is clearly distinguished. Comparing the different angular resolution 
images, the flux and size of MM1 and MM2 decrease when angular resolution increases. 
In the images with robust <0. 5 (synthesized beam size of about 1"), MM1 is unresolved 
and MM2 appears elongated in an E-W direction joining MM2(E) and MM2(W). A two- 
Gaussian fit on MM2 in the images with robust<0.25 shows MM2(E) marginally resolved 
(< 0.8") with a P.A. near to 50° ± 15°, while in the robust=0 image it appears unresolved. 
On the other hand, MM2(W) appears spatially resolved in both images (robust=0.25 and 
robust=0 images), with sizes and orientations similar to those observed in the very extended 
configuration image (see Table [3]). The two components of MM2 are only well resolved in 
the highest angular resolution image (image with robust=0). 

The fluxes of MM1 and MM2 decrease in a different way with the angular resolution 
of the images. We note here that the area of the synthesized beam changes approximately 
a factor of 4 between the images with robust=2 and robust=0.5, and a factor of 8 between 
the images with robust=0.5 and robust=0 (see the caption of Fig. [3]). The flux of MM1 
decreases by 22% between the robust=2 and the robust=0.5 images, and by only 6% between 
the robust=0.5 and the robust=0 images. Thus, the MM1 flux decreases by 28% (170 mJy) 
in all, with respect to the flux obtained with the robust=2 image. In the case of MM2, the 
flux decreases by 37% between the images with robust 2 and 0.5, by 27% between those with 
robust 0.5 and 0.25, and by 12% between those with robust 0.25 and 0. Considering the 
flux of MM2(E) and MM2( W) altogether (each one contributes about half of the total flux 
measured in MM2), the total decrease of the MM2 flux is of about 76% (340 mJy). Finally, 
these images also show about 135 mJy of extended emission, which is seen as weak emission 
bridges between MM1 and MM2 (see for instance the image with robust=0.5, Fig. [3]), and 
as a weak filamentary structure towards the west of MM2 (Fig. [3j image with robust=0.25). 

The flux variations with resolution observed in MM1 and MM2 indicate the existence of 
extended emission associated with these sources. These results also indicate that the emission 
of MM1 comes mainly from a compact component embedded in a weak extended envelope. In 
the case of MM2, most of the emission is extended and elongated in an East- West direction; 
the extended envelope surrounds the compact MM2(E), and the resolved source MM2(W). 
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In section §5 we discuss the morphology and the flux variations of the sources with the 
angular resolution as a possible evidence of MM1 and MM2 being in different evolutionary 
stages. 

GGD27 has been previously observed at 1.4 mm with the BIMA and at 1.33 mm with 
the SMA interferometers (Gomez et al. 2003; Qiu & Zhang 2009) and at 450 /iin with 
the JCMT (Jenness, Scott, & Padman 1995; McCutcheon et al. 1995; Thompson et al. 
2006). From BIMA maps (using a similar angular resolution as our compact configuration 
image), Gomez et al. (2003) reported a total flux of 1.18 Jy, in good agreement with the 
flux measured here in the compact configuration image. Qiu & Zhang (2009) reported a 
total flux of 2 Jy measured with the SMA at 1.33 mm, which is a factor of 2 greater than 
that reported here. The data used by Qiu & Zhang (2009) included SMA observations with 
compact and subcompact configurations, and obtained a better uv-coverage and a greater 
sensitivity to extended structures. The different JCMT observations at 450 fim showed a 
dusty core of ~ 40" of diameter nearly centered on MM1 and MM2. The flux obtained by the 
JCMT observations is of 97±30 Jy (Thompson et al. 2006), while the total flux recovered in 
our SMA 456 fim image is 24.3 Jy, which is 25% of the flux measured with the JCMT. Since 
the SMA in its compact configuration is sensitive to structures smaller than 22" at 1.36 mm 
and 7" at 456 fim, the detected sources at both wavelengths correspond to the emission of 
compact structures located inside the core detected with the JCMT. 

The SMA observations that we present here show that at 1.36 mm the continuum flux 
of MM1 is higher than that obtained for MM2, while at 456 /im, MM2 is the strongest source 
in the region. The spectral indices of MM1 and MM2, using a power-law model between 
1.36 mm and 456 fim, are 2.3±0.5 and 3.6±0.4, respectively. We find that the free-free 
contribution based on 3.6 cm emission extrapolated to 1.36 mm is negligible (<5 mJy for 
MM1 and <0.5 mJy for MM2 at 1.36 mm, assuming free-free spectral indices of 0.2 and 0.6 
respectively; see §3.2 below). Up to now, MM1 has been considered as the main source in the 
region and thus the one that contributes most of the bolometric luminosity, obtained from 
IRAS and low angular resolution (mostly single-dish) observations (see Fig. [5]). However, 
the 456 fim SMA observations suggest that MM2 contributes substantially to the total 
luminosity of the region. 



3.2. Spectral Energy Distributions 

The high angular resolution data obtained with the SMA towards GGD27, presented 
in this paper, allow the extension of the Spectral Energy Distributon (SED) of the sources 
MM1 and MM2. MM1 has been observed with high angular resolution in a very wide range 



13 



of wavelengths (from radio to NIR; see Table H|). MM2, on the other hand, has only been 
detected at three different wavelengths: the two reported here and at 3.5 cm by Gomez et al. 
(1995). In this section we analyze the nature of both sources by fitting the SED of MM1 
and deriving the spectral index of MM2 in the (sub) millimeter range. 

The SED of MM1 was built using the data from radio (VLA) and (sub)millimeter (BIMA 
and SMA) interferometric observations, and from infrared observations with good angular 
resolution (see Table H]). The selected data have an angular resolution comparable to our 
456 /im observations (i.e., about 2"). In this way, the measured emission comes from the 
vicinity of MM1, avoiding the contamination due to emission from other sources. Figure 
H] shows the SED of MM1 and a fit to these points. The fit minimizes the \ 2 following 
the algorithms described in Canto, Curiel, & Martinez-Gomez (2009). The function fitted 
consists of a modified black body plus free-free emission model: 

S v = B v (T dust )(l - e- T »)n s + a v v a , 

where B v is the Plank function with a temperature T^ust, &s is the solid angle subtended 
by MM1, t u is the opacity which depends on frequency as r u = r • (u/uq) , ft being the 
dust opacity index and Uq an arbitrary reference frequency, and a v = cio/^o- Since this 
model includes a single modified black body, the points with A < 10 /im shown in Figure 
H] were not taken into account in the fit. We could fit these NIR data with an additional 
modified black body, but the uncertainty on the ft parameter prevents us from using a second 
component. However, these NIR data points seem to be associated with dust emission at 
high temperatures. Table [5] shows the results of the best fit with 6-independent parameters 
{Tdust, T o, 0i a o an d a), as well as the uncertainties provided by the algorithm. The 
obtained spectral index of the free-free emission (a = 0.18 ± 0.01) is consistent with the 
value of a = 0.2 obtained by Marti et al. (1993) and by Gomez et al. (2003) independently, 
implying a negligible free-free contribution in the (sub)millimeter range. The obtained dust 
opacity index (ft = 0.53 ± 0.05), agrees well with the value reported by Gomez et al. (2003) 
(ft = 0.6), obtained from a power-law fit from the flux densities obtained at 7, 3.5 and 
1.4 mm. The low value of the dust opacity index suggests the presence of processed dust 
grains in MM1 (Draine 2006 and reference therein). In addition, the dust temperature 
obtained for MM1 is 109 ±3 K and the estimated size of the source is 1.1 ± 0.4 arcsec 2 , 
consistent with the compact nature of this source. Finally, the bolometric luminosity of 
MM1 was estimated as 3300 L . This luminosity is smaller than the luminosity obtained 
from low angular resolution data (~2x 10 4 L ; see Fig. [5]), and it could be considered as a 
lower limit, since the high angular resolution observations used in this fit do not include the 
outer parts of the main core (only includes the dense structure surrounding MM1), which 
are also heated by this source (see also the discussion in section §4). 
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In the case of MM2, there are only three high angular resolution (< 1") flux estimates 
(Fig. [5]): the fluxes at 1.36 mm and 456 /zm (this paper), and the flux at 3.5 cm (Gomez 
et al. 1995). The spectral index between 1.36 mm and 456 /im (a = 3.60 ±0.4) corresponds 
to a dust opacity index f3 between 1.2 and 2.0 (accounting for the error bars) which is similar 
to the averaged opacity index found in the interstellar medium ((3 = 1.6; Draine 2006). This 
could indicate that the dust in MM2 contains mainly small particles, such as those found 
in the interstellar medium, which suggests that MM2 could be a very young source. In 
addition, the 3.5 cm emission is consistent with free-free emission, since the expected dust 
emission at this wavelength (using a spectral index of 3.6) is only 6% of the measured radio 
continuum flux. This free-free emission could be arising from an ionized outflow associated 
with a young intermediate/high-mass protostar. 



3.3. Masses and Column Densities 

The equations used to estimate the mass and the H 2 column density from the continuum 
dust emission, expressed in appropriate units are (Hildebrand 1983): 

/ 0.048^- -. \ t--2 q 

M[M ] = 3.2643 • 10 6 • R ■ [C / ,~ * W 

iV^[cm- 2 ] = 8.4208 -10 23 -^y (2) 

where R, T d , D, S u , k u and fl are the gas-to-dust ratio (which is assumed to be 100), the 
dust temperature in Kelvin, the distance to the source in kpc (adopted as 1.7 kpc), the 
measured flux density in Jy, the dust opacity in cm 2 • g and the angular size in arcsec 2 . 
The frequency v is in GHz. In order to obtain the dust opacity, h\^ mm = 1.0 cm 2 -g _1 was 
assumed for thin ice mantles (Ossenkopf & Henning 1994), scaling this value with frequency 
as z/ 3 , where (3 is the dust opacity index derived in §3.2 for MM1 and MM2. The hypothesis 
of thin ice mantles towards GGD27 is suggested by the detection of absorption features 
of water ices at 3 and 6 /im, and CO2 ices at 4.3 and 15.2 /im (Peeters et al. 2002; van 
Diedenhoven et al. 2004). The mass, column density and number density were obtained for 
each source (see Table E]), from the measured flux densities at 1.36 mm (high and low angular 
resolution images) and at 456 /im. For MM1, we use T<2 US i=109 K and /3=0.53 derived from 
the SED fit in §3.2, while for MM2 we adopt (3=1.6 (see §3.2) and Tdust=35 K as derived 
from observations of several CH3CN (12-11) transitions associated with this source (Qiu & 
Zhang 2009). The uncertainties were estimated from the errors in Tdust, P and the flux 
densities. These uncertainties do not take into account the effects introduced by the opacity, 
which could increase the value of the derived masses. 
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Using the flux densities obtained with the two-Gaussian fit in the low angular reso- 
lution image at 1.36 mm, the masses estimated for MM1 and MM2 are 5.6±0.2 M and 
14.9±0.4 M Q , respectively. As expected, these values are similar to those obtained using 
the 456 /im flux densities (3±1 M for MM1 and 17±2 M for MM2). On the other hand, 
from the high angular resolution data at 1.36 mm, we estimate that the mass of MM1 is 
4.1±0.2 M and the masses of MM2(E) and MM2(W) are 1.9±0.2 M Q and 1.7±0.3 M , 
respectively. Thus, the extended dust emission surrounding MM1 and MM2 contains about 
14 M , most of which (~ 12 M ) seems to be associated with MM2. The protostars in MM2 
still have a huge reservoir of gas and dust to accrete, while the gas and dust supply for the 
protostar in MM1 is almost an order of magnitude smaller (see discussion in section §4). 

The total mass obtained from the fluxes in the low angular resolution image (~20 M ) is, 
approximately, 5% of the total mass obtained for the ~arcmin dusty core detected by single- 
dish observations (e.g., Hunter et al. 2000), within which MM1 and MM2 are embedded. 
Thus, only a small fraction of the material in the large scale core is associated with the 
compact structures. This suggests that the main core still has enough material to form a 
cluster of intermediate/high-mass protostars or to contribute with more accretion into MM1 
and MM2. 

Finally, based on the high angular resolution observations at 1.36 mm, we obtain a 
lower limit for the column density of MM1 (1.13 x 10 25 cm" 2 ), MM2(E) (5.5 x 10 24 cm" 2 ) 
and MM2(W) (8.9 x 10 23 cm -2 ). In addition, we also estimate lower limits for the volume 
number densities: 1.85 x 10 9 cm" 3 for MM1, 9.0 x 10 8 cm" 3 for MM2(E) and 1.1 x 10 8 cm" 3 
for MM2(W). Hence, the volume number density towards MM1 and MM2(E) could be one 
order of magnitude greater than in MM2(W). Such high densities (> 10 8 cm -3 ) have also 
been found in other high-mass star formation regions (e.g. IRAS20126, Cesaroni et al. 2005; 
Cep A, Jimenez-Serra et al. 2007; W51 North, Zapata et al. 2009), around circumstellar 
structures (possible accretion disks). 



4. DISCUSSION 

The high angular resolution observations at 1.36 mm towards GGD27 show that MM1 is 
a compact source, which spatially coincides with the powerful radio continuum jet observed 
in this region (Marti et al. 1993, 1995; Gomez et al. 1995). MM1 appears unresolved after 
applying a uniform weighting to the very extended configuration data at 1.36 mm. The SNR 
of MM1 in this image is ~ 100 (see image with robust in Fig. [3]), so taking into account 
some additional considerations, it is possible to obtain a more restrictive upper limit for the 
diameter of this source than the / .'36 (i.e., about 600 AU at 1.7 kpc) given in §3.1. To make 
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this new estimation, the dust emission from MM1 is assumed to arise from a flat envelope 
in the shape of an accretion disk. This hypothesis is based on the presence of a radio jet 
and a molecular outflow asociated with MM1. Such a disk should be oriented perpendicular 
to the jet's axis, whose P.A. is 21° (Marti et al. 1993). On the other hand, the uniform- 
weighted image obtained with the extended configuration data has a synthesized beam of 
0'/72 x 0'/31, with the major axis of the synthesized beam (P.A.=19?4) oriented almost in the 
same direction as the radio jet. Thus, the major axis of the synthesized beam follows the 
direction of the radio jet, while the minor axis of the synthesized beam follows the orientation 
of the putative accretion disk. Therefore, given the high SNR in MM1 and assuming that 
the compact dust emission comes from an accretion disk, a reasonable upper limit for the 
diameter of this disk would be half the minor axis of the synthesized beam, which is 0'.'16 or 
272 AU at the assumed distance. If the disk hypothesis is right, this is one of the smallest 
disks found around a high-mass protostar. Some of the disks detected around protostars 
with 10 3 — 10 5 L have a radius ranging between 70-1000 AU and a mass between 1-20 M Q 
(Cesaroni et al. (2007); Patel et al. (2005); Shepherd, Claussen, & Kurtz (2001); Rodriguez, 
Zapata, & Ho (2007); Franco-Hernandez et al. (2009); Galvan-Madrid et al. (2010).). 

In addition, although the diameter of the disk in MM1 is similar (perhaps somewhat 
greater) to that of disks around low-mass protostars, its mass (~4 M ) is about two orders 
of magnitude larger than those found in low-mass protostars (~ 0.03 M ; e.g, Palau et al. 
2006). These massive disks have been reported in other sources. The existence of a disk 
around MM1 would imply a star formation mechanism based on the accretion of circumstellar 
material. 

The morphology of the northern source (MM2) varies with the angular resolution of the 
1.36 mm images (Fig. [3]). These images show that MM2 contains a compact core, MM2(E), 
and an extended structure, MM2(W), both surrounded by an extended and massive envelope. 
MM2(E) is probably a compact core with a central protostar as indicated by the presence 
of a weak radio continuum source and a water maser (Gomez et al. 1995; Marti et al. 1999; 
Kurtz & Hofner 2005). MM2(W) is probably a prestellar core with enough mass to form 
an intermediate-high mass protostar (probably an early-B type). Spectral line observations 
with high angular resolution may show whether MM2(W) is collapsing or not. In addition, 
given the small projected distance between MM2(E) and MM2(W) (about 1600 AU), the 
two sources might be gravitationally bound. There is enough mass to allow orbital motions 
of < 1.5 km s -1 . Finally, MM2 is characterized by strong submillimeter emission (19 Jy at 
456 /mi), low temperature (T ~ 35 K, Qiu & Zhang 2009), high density (> 10 7 cm" 3 ), a dust 
opacity index close to that of the interstellar medium (j3 ~ 1.6), non-detection in the MIR 
and NIR, association with a young molecular outflow (Qiu & Zhang 2009) and a massive 
dusty envelope (11 M Q ). These properties are similar to the characteristics of the Class 
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low-mass protostars (Andre, Ward-Thompson, & Barsony 1993; Andre & Montmerle 1994). 
However, MM2 seems to be embedded in an envelope with a mass four times greater than 
the envelopes of Class low-mass protostars (which have 1-4 M Q ; e.g., NGC 1333 IRAS4A, 
Girart, Rao, & Marrone 2006; L723, Girart, Rao, & Estalella 2009; IRAS16293A, Rao et al. 
2009). The higher mass of the envelope could also be explained by multiple protostars within 
it. 

MM1 and MM2 show different characteristics. MM2 has a (3 ~ 1.6 similar to the 
interstellar medium (Draine 2006 and references therein), while MM1 has a /3 ~ 0.5 which 
suggests more processed dust (Beckwith, Henning, & Nakagawa 2000; Draine 2006; Natta 
et al. 2007). The temperature of MM1 is greater than that of MM2 and the free- free 
emission of MM1 is two orders of magnitude greater than that of MM2. In addition, the 
outflow associated with MM1 has a mass of 407-460 M Q and a dynamic time scale of ~ 10 5 yr 
(Yamashita et al. 1989; Ridge & Moore 2001), while the outflow associated with MM2 has a 
mass of 0.22 M and a dynamic time scale of ~ 2 x 10 3 yr (Qiu & Zhang 2009). Finally, the 
mass of the envelope of MM1 is almost an order of magnitude smaller than that of MM2, 
which could indicate that the protostars in MM2 are at an earlier stage of accretion from 
the envelope as compared to MM1. All these characteristics are consistent with MM2 being 
in a younger evolutionary stage than MM1. 

MM1 is the strongest source at 1.36 mm, while MM2 dominates the emission at 456 /im. 
Therefore the present SMA data indicate that both sources significantly contribute to the 
total luminosity of the region, unlike previous publications (e.g., Rodriguez & Reipurth 
1989; Marti et al. 1999; Molinari, Noriega-Crespo, & Spinoglio 2001; Qiu et al. 2008), which 
pointed at MM1 as the source responsible for all the bolometric luminosity in the region. 
With the available data, it is unclear which is the more luminous source, but their SEDs 
(Fig. E]) suggest comparable contribution to the total bolometric luminosity (estimated as 
~ 2 x 10 4 L from IRAS, JCMT and CSO data; see Fig. EJ. According to Table EJof Molinari 
et al. (1998), a luminosity of 10 4 L Q would suggest a 13-18 M Q protostar. 

The luminosity of MM1 estimated from the fit to its SED is 3300 L , which is smaller 
than half the bolometric luminosity of the region. However, it has been shown (e.g., Yorke & 
Bodenheimer 1999) that the inclination of a disk-protostar system is an important parameter 
for the estimations of the luminosity from the SED. In this case, the radio jet associated with 
MM1 is almost in the plane of the sky (Aspin et al. 1991; Marti et al. 1995), which implies 
that the circumstellar disk-like structure is approximately edge-on. The large column density 
derived in the disk-like structure yields a visual extinction of at least 1.4 x 10 4 . This implies 
that along the line of sight most of the mid-IR radiation is effectively absorbed by the disk. 
In terms of total luminosity, this suggests that a significant fraction of the radiation escapes 
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through the outflow cavity (Yorke & Bodenheimer 1999; Krumholz et al. 2009). Therefore, 
the bolometric luminosity of the young star in MM1 is likely larger than the measured 
3300 L , and it would likely correspond to a Bl ZAMS star. However, the presence of the 
massive disk and of the powerful outflow indicates that it is still in a very active accretion 
phase, so the protostar could become a BO ZAMS star, or even a more massive star. 

Finally, assuming that MM1 is excited by a Bl-type protostar, it should exhibit a 
compact HII region. Free-free emission at several wavelengths has been reported towards this 
source, but the spectral index (0.18±0.01, section §3.2) clearly corresponds to that expected 
for an ionized wind. This suggests that the protostar in MM1 has not yet developed an HII 
region and therefore it could be considered as a true high-mass protostar. 

5. CONCLUSIONS 

The present paper presents a study of the GGD27 high-mass star formation region based 
on SMA observations of the continuum emission at 1.36 mm (compact and very extended 
configurations) and 456 /im (compact configuration). These observations show compact 
dust emission in two regions, MM1 and MM2. Using these new observations together with 
previously published results, we analyze the SEDs of MM1 and MM2. The masses and 
densities of these two sources are also estimated. The main conclusions of this paper are: 

1. The dust continuum emission at 1.36 mm and 456 /im comes from two compact 
structures (MM1 and MM2), separated 7" (about 12000 AU in projected distance) and 
embedded in a ~ 40" core, that was detected by single-dish observations. The high angular 
resolution observations at 1.36 mm show that MM1 is unresolved (with an estimated radius 
of less than 150 AU), with an estimated mass of 4.1 M . MM1 spatially coincides with 
an IRAS source of 2 x 10 4 L and it is associated with a powerful molecular outflow and 
a radio continuum jet. These results suggest that the dust emission in MM1 probably 
arises from a massive accretion disk. At subarcsecond angular scales, MM2 splits into two 
main components, MM2(E) and MM2(W), that are separated by 1" (1700 AU in projected 
distance), each with an estimated mass of ~ 2 M Q . MM2(E) is probably a compact core 
with a young protostar inside (indicated by a weak radio continuum source and a water 
maser) and MM2(W) is probably a prestellar core. Both sources are embedded in a massive 
envelope of at least 11 M . 

2. The characteristics of MM1 and MM2 are very different: (a) MM2 has a steeper 
spectral index at (sub) mm wavelengths, suggesting a less processed dust than in MM1, (b) 
MM1 has a higher temperature, suggesting that there is a higher mass protostar inside MM1 
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than in MM2, (c) the molecular outflow associated with MM2 is younger and less massive 
than the one associated with MM1, (d) MM1 is unresolved even at subarcsecond scales, while 
MM2 appears extended with two main components and faint emission surrounding them. 
These characteristics suggest that MM1 and MM2 are in different evolutionary stages, with 
MM2 younger than MM1. 

3. The SMA observations show that both MM1 and MM2 contribute significantly to the 
bolometric luminosity of the whole region (~ 2x 10 4 L Q ). A fit to the SED of MM1 indicates 
that a lower limit to the luminosity of this source is 3300 L , which implies that MM1 is 
at least a Bl ZAMS protostar. However, the likely continued accretion in MM1 indicates 
that this source will accrete several solar masses. At present it is not clear what kind of 
protostar is inside MM2. However, these SMA observations suggest that MM2 contributes 
an important fraction of the bolometric luminosity of this region. 
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A. APPENDIX 

The phase transfer calibration method relies on the fact that phase variations at high 
frequencies may be correlated in some manner to those at low frequencies, where the SMA 
is more stable and more sensitive (the phase noise from the system and the sky are much 
better when working at such low frequencies). Since the atmosphere is almost non- dispersive 
(i.e., its refraction index does not depend strongly on frequency, away from resonant lines), 
to first order, the ratio of phases at two different frequencies is expected to be linear with the 
frequency ratio (i.e., oc see e.g. Hunter et al. 2005). Phase transfer is expected to be 

useful with sub millimeter interferometers as the SMA, due to the lack of strong calibrators at 
high frequencies, which could make the standard techniques of phase calibration unavailable. 

VX Sgr, the evolved star included in the observing cycle, has maser emission from SiO at 
215.596 GHz and from H 2 at 658.007 GHz. Both masers were detected simultaneously by 
the SMA, so we analysed their phase behavior in order to check the phase transfer procedure. 
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In first place, the beginning and the end of the track were flagged due to poor quality (when 
elevation of GGD27 and VX Sgr were below 30°). The data were corrected for variations in 
system temperature and bandpass, and as in §2.2.2 we made a pseudo-continuum channel for 
each maser. The phases of the masers showed in general good linear correlations (above an 
80%). These correlations gave different values for the slope and the abscisa offset depending 
on the baseline lengths (panels (a), (b) and (e) of Fig. [6]). The linear dependence of the 
slopes of the correlation with the length of each baseline suggests that there are instrumental 
effects between the antennas and the correlator. On the other hand, the phase drifts seem to 
be almost linear during long periods of time. Hence, after self calibrate each set of maser data 
with a time averaging of 10 minutes, we removed these slow drifts. After that we tried to 
cross-correlate the phases of both masers again, but now they were uncorrelated. The panels 
(c) and (d) of Figure [6] show the behaviour of the phases of both masers for the two baselines. 
Ideally, we expect the phases to be linearly related by the ratio between frequencies, but in 
this case, the rms of the phase at 658 GHz (about 40°), is probably masking the theoretical 
behaviour. This would occur at high frequencies because of the higher rms of the phases. 

Despite the high rms of the phase at high frequency we tested the phase transfer tech- 
nique. We scaled up the 215 GHz phase gain solutions to that of the 658 GHz uncalibrated 
data of VX Sgr, using the Miriad task PHATRANS. PHATRANS extrapolated the phase 
corrections from the phase gains determined at 215 GHz with the SiO maser, to the phase 
corrections at 658 GHz of the VX Sgr uncalibrated data, using the linear relation between 
the phase gains of the SiO and water pseudo-continuum channels at both frequencies. Then, 
we performed the flux calibration using Ceres. After continuum subtraction, the output 
cleaned and restored velocity cube image of the water maser at 658 GHz (spectral resolu- 
tion of 0.37 km s -1 ) had an rms of 18 Jy beam -1 , and a SNR of 22. Five iterations of self 
calibration on this velocity cube, improved its SNR. The final velocity cube has an rms of 
16 Jy beam -1 and a SNR of 36. 

To check the validity of the phase transfer calibration, we also performed a similar cal- 
ibration of the data of VX Sgr at 658 GHz applying directly the phase gain solutions of the 
the pseudo-continuum channel of the water maser. We refer to this calibration technique 
as phase referencing. The only difference with the phase transfer procedure was the phase 
calibration step. In this case, the final velocity cube image of the maser had an rms of 
16 Jy beam -1 and a SNR of 38. Self calibration did not improve this image substantially 
(less than 0.5 in SNR). The phase transfer calibration of the water maser of VX Sgr shows a 
similar output data than the phase referencing calibration: a) The SNR of the velocity cubes 
obtained with phase transfer and phase referencing are similar (36 and 38, respectively), b) 
Both velocity cubes have the same qualitative information (i.e., the same channels show 
emission over 3 a). c) The spectra obtained with the two calibration techniques have qual- 
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itatively the similar shape (same width and peak position). The difference between both 
spectra is smaller than 3 o for practically all the spectral channels, d) The moment images 
are similar qualitatively, showing an unresolved source at the phase center, e) The total flux 
measured in the moment images is only 5% smaller in the phase transfer image. The SNR 
of the phase tranfer and the phase referencing images are 19 and 20, respectively. Thus, the 
analysis indicates that in this case, phase transfer could give a similar output data (images 
and spectrum) than the phase referencing, taking into account the rms of the images and 
the calibration uncertainties. 

We also tried to apply the phase transfer procedure to the GGD 27 data. However, in 
this case, the source was undetected, preventing to start with self calibration. Therefore, our 
attempt to apply the phase transfer technique seems to be quite successful when calibrating 
the bright water maser of VX Sgr, but it failed with a fainter source such as the continuum 
emission of GGD 27. We can conclude that the phase transfer technique described here 
works, at least in this case, for a source strong enough to apply further self calibration. 
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Fig. 1. — 1.36 mm (220 GHz) continuum images of the region. The left panel shows the 
low angular resolution image and the right panel shows the high angular resolution image, 
both obtained with natural weighting. Contours are -5, -3, 3, 5, 10, 20 and 40 xlO.O mJy 
beam -1 , the rms noise of the low angular resolution data and -5, -3, 3, 5, 10, 30, 75 and 
120 x3.0 mJy beam -1 , the rms noise of the high angular resolution data, respectively. The 
filled squares mark the positions of the sources detected at 3.5 cm (Gomez et al. 1995), the 
crosses mark two water masers (Gomez et al. 1995; Kurtz & Hofner 2005) and the triangle 
marks the position of a CH3OH class I maser (Kurtz et al. 2004). The positional accuracy 
of these sources is O'/l, except for the water maser located to the north of MMl, on the radio 
jet tip at 3.5 cm (see Fig. [2]), which is O'.'Ol (Kurtz & Hofner 2005). The synthesized beam 
are shown as a gray solid ellipse at the bottom left corner of each image. The dashed circle 
shows the full width half maximum contour of the primary beam (which has a radius of 25"). 
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Fig. 2. — Overlay of our 456 /mi (658 GHz) continuum image (black contours) with the 
3.5 cm continuum VLA image (red contours) obtained by Gomez et al. (1995) towards the 
central region of GGD27. Contour levels are -6, -3, 3, 6, 9, 12, 15, 20, 30, 40 x457 mJy 
beam -1 , the rms noise of the 456 /im image and -5, 5, 6, 7, 8, 10, 12, 15, 20, 40,80, 120, 
160, 200 x9.1 /iJy beam -1 , the rms noise of the 3.5 cm image (the latter, as appeared in 
the original article). The 456 /im image is presented without primary beam correction. The 
synthesized beams are shown at the bottom right corner and the primary beam is indicated 
by a dashed circle with radius 8'/5 centered on the phase center. Symbols are the same as 
in Fig. [TJ The dust emission towards MM1 coincides with the radio continuum jet, the 
suspected powering source of the HH80-81 bipolar system. Likewise, the dust emission in 
MM2 coincides with a very weak radio continuum source (the smallest contour around the 
cross shows the 5 a emission at 3.5 cm). The weak source to the north, near the primary 
beam edge is an artifact of the cleaning process. 
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Fig. 3. — Combined images of the extended and compact configuration at 1.36 mm (220 GHz) 
obtained with different values of the robust parameter. The value of the robust weighting is 
indicated within each panel. Contours are -5, -3, 3, 5, 10, 20 and 40 x8.1, x5.3 mJy beam" 1 , 
the rms noise of the images with robust 2 and 0.5 (top panels), and -6, -3, 3, 6, 12, 30 and 
60 x4.2 and x3.9 mJy beam -1 , the rms noise of the images with robust 0.25 and (bottom 
panels), respectively. Symbols are the same as in Fig. [TJ Synthesized beams are 5'.'0 x 2'.'4 
(36?3), 2"0 x 1*3 (45?0), 0"9 x 0"7 (31?0) and 0'.'8 x 0'/4 (15?1) for the images with robust 2, 
0.5, 0.25 and 0, respectively, and are shown as gray solid ellipses at the bottom left corner 
of each image. 



- 30 - 



Wavelength (mm) 
100 10 1 0.1 0.01 




1 10 100 1000 10000 100000 

Frequency (GHz) 

Fig. 4. — Spectral Energy Distribution from radio to NIR wavelengths for MM1 (filled circles with error 
bars). The point-dashed curve is the free-free emission and the blue curve the modified black body fit of 
MMl's SED (see text). The three NIR points taken from Table [4] with wavelengths below 10 /xm, are not 
included in the fit due to the uncertainty on the /3 parameter of an additional modified black body. 



Table 1. Continuum observations 



Continuum 


Frequency 


Date 




Synthesized beam 
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rms 


Map 
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P.A. 
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(") 


(°) 


(") 


(mJy beam -1 ) 


Lower Side Band 
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24 
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34.2 


50 


11 


UpperSide Band 
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10 
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8.1 x 3.0 


34.1 


50 


10 
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2007 May 29 


0.70 x 0.39 


13.2 


50 


4.3 


UpperSide Band 
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2007 May 29 


0.67 x 0.36 
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LSB + USB 


220 
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0.68 x 0.38 


13.0 


50 


3.0 
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24 
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( Q ) Primary beam corrected image. 
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Fig. 5. — Spectral Energy Distribution from radio to NIR wavelengths for MM1 (red circles) and for MM2 
(blue circles). The red solid curve accounts for the fit of MMl's SED (see text and Fig. [4J. The blue dashed 
curve is a power law adjusted to the mm and submm points of MM2 with spectral index 3.6, corresponding 
to a (3 of 1.6 (see section §3.2). The green circles correspond to the low resolution data taken from literature 
(single-dish and IRAS measurements). Most of the error bars are smaller than the size of each symbol. 
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Fig. 6. — (a) and (b): 215 GHz (red circles) and 658 GHz (blue circles) uncalibrated phase of the masers of 
VX Sgr from two baselines with very different separations (the 3-4 baseline with 14.5 m and the 5-6 baseline 
with 64.0 m). Each point represents the average value of a 30 seconds integration, (c) and (d): 5 minutes of 
the phase behaviour at 215 GHz (red circles) and 658 GHz (blue circles) for both baselines, after subtracting 
a first order polynomial fitted to the raw data. Error bars are derived from the rms noise of this 5 minutes 
data, (e): cross correlation slope versus the baseline separation. The straight line is a fit to most of the 
measurements (filled circles) but the correlations of baselines 2-3, 3-5 and 4-6 (which are located well below 
the general trend) . 
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Table 2. Flux calibration. Planets, moons and other calibrators. 



Source 


Frequency 


Measured Flux< a > Expected Flux' 6 ) 


Deviation 


Offset^) 


Observing Date 




(GHz) 


(Jy) 


(Jy) 


(%) 


(") 




j ly 1 i-zui 


01 

Z10 


2.13 


(d) 




n no/i n aii 


AUg Z4,ZUU0 


j ly 1 i-zui 


99£ 


2.18 


2.37 


o.U 


n non n nna 
U.U^U,-U.UUo 


/\Ug Z4,ZUU0 


'iCAKA Q 
OU4D4.0 




12.87 


(d) 




U.U1U,U.U40 


/\Ug Z4,ZUUD 




225 


28.27 


27.85 


1 ^ 


u.uuo,u.uou 


Alio- 94 900^ 


Ganymede 


215 


3.06 


3.52 


12.9 


0.022,0.015 


Aug 24,2005 


Ganymede 


225 


3.29 


3.84 


14.3 


0.010,0.013 


Aug 24,2005 


CallistoM 


215 


3.34 


3.57 


6.2 


-0.003,-0.010 


Aug 24,2005 


Callisto( e ) 


225 


3.63 


3.90 


6.7 


0.005,0.014 


Aug 24,2005 


Neptune 


215 


13.80 


13.41 


3.0 


-0.019,-0.001 


Aug 24,2005 


Neptune 


225 


14.21 


13.93 


2.4 


0.014,0.003 


Aug 24,2005 


Ceres 


215 


0.83 


0.99 


15.9 


0.030,0.016 


Aug 24,2005 


Ceres 


225 


0.83 


1.08 


23.0 


0.005,0.012 


Aug 24,2005 


J1911-201 


658 


... (/) 


< 1.5 






Aug 24,2005 


3C454.3 


658 


13.50 


... (d) 




-0.004,0.013 


Aug 24,2005 


Ganymede 


658 


25.82 


29.40 


12.4 


-0.057,-0.057 


Aug 24,2005 


Callisto 


658 


32.38 


30.38 


6.4 


0.028,0.054 


Aug 24,2005 


Neptune 


658 


46.30 


68.84 


32.8 


-0.008,0.056 


Aug 24,2005 


Ceres( £ ) 


658 


9.03 


8.67 


3.9 


0.032,0.009 


Aug 24,2005 


J1911-20l( e ' 


215 


1.37 


1.24 


10.5 


0.004,-0.004 


May 29, 2007 


J1911-20l( e ' 


225 


1.29 


1.16 


11.2 


0.004,-0.002 


May 29, 2007 


3C454.3< 9 ) 


220 


7.69 


7.41 


3.8 


0.000,0.002 


May 29, 2007 


mwc349 


215 


1.44 


... (d) 




0.002,-0.003 


May 29, 2007 


mwc349 


225 


1.92 


1.73 


11.0 


-0.000,-0.000 


May 29, 2007 



( a )Flux density within a box surrounding the emission from the calibrator. 

( & )Data interpolated from the SMA monitoring program of calibrators and from the planetary visibility function 
calculator of the SMA webpage (Gurwell, M.). 

( c ) Peak position from a fit performed with the MAXFIT task of Miriad. The offsets are relative to the phase 
center of each source. 

( d ) Data not available from the SMA monitoring program of calibrators. 
( e 'Flux calibrator. 

(■^Undetected at this frequency. 

(s^The observed data at 220 GHz is an average from data at 215 and 225 GHz (7.69 and 9.27 Jy respectively). 
This value is compared with that measured by the SMA monitoring program at 219.8 GHz on the same day of our 
observations. 
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Table 3. Source parameters 



Frequency 
(GHz) 


Beam 


Source 


lgh 19 m 


DEC( a ' 
-20° 47' 


(mjy beam -1 ) 


(mjy) 


Size' 6 ) 
(") 


P.A. 

(°) 


220 


8'.'1 x 3'.'0; 34? 1 


MM1< C ) 


12?114 


30'.'61 


455±4 


610±8 


3.5 ±0.1 x 1.3 ±0.2 


176±3 


220 


8'.'1 x 3'.'0; 34? 1 


MM2< C ) 


12:488 


27'.' 30 


260±4 


443±9 


4.9 ±0.2 x 1.5 ±0.5 


164±3 


658 


2'.' 4 x l'.'l; 37? 1 


MM1 


12! 105 


30'.' 90 


4900±900 


5300±1700 


$ 1-17 




658 


2('4 x l'.'l; 37? 1 


MM2 


12?517 


27'.'09 


18600±900 


19000±1500 


$ 1-17 




220 


0'.'68 x 0"38; 13?0 


MMlW 


12?099 


30" 77 


409±3 


441±6 


< 0.16( e > 




220 


0'.'68 x 0{'38; 13?0 


MM2(E)( d ) 


12?515 


27'.' 23 


47±3 


57±6 


<0.36 




220 


0'.'68 x 0"38; 13?0 


MM2(W)( d > 


12:456 


27'.'00 


17±3 


50±10 


1.1 ±0.2 x 0.5 ±0.1 


34±11 


Note. — Natural weighting was applied in all the images. Fits were carried out with the IMFIT routine of AIPS. The errors are obtained 
from the fits and do not include the 20% of the calibration uncertainties. The fluxes from the image at 456/^m are primary beam corrected. 



( a )j2000. 



O Deconvolved sizes. When deconvolution fails wc use half the synthesized beam size area ((1/2) ■ (Trb ma j ■ b m i n /4)) as an upper limit for 
the size of the region and calculate its diameter as ■ Beam arca/-7r. This could be a useful effective diameter for elliptical synthesized 
beams with high eccentricity. 

( c 'ln the low angular resolution image, the fit was carried out with two Gaussians simultaneously, because MM1 and MM2 share part 
of their emission. The uncertainties on the flux density and the sizes could be higher than the given by the fitting algorithm because the 
two sources are not well resolved. 

( d ) ln the high angular resolution image, the fit was carried out with one Gaussian for MM1 and independently, two Gaussians for MM2. 

( e ) Scc §4 for the explanation of this upper limit. 
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Table 4. MM1 SED data 



Wavelength 


Integrated Flux 


Error 


Aperture 


Reference 


mm 


Jy 


Jy 


// 




60 


0.00272 


0.00002 


0.8 x 0.5 


1 


36 


00296 


00002 

V 7 • V 7 V 7 V 7 V 7 


5x03 


1 


35 


0.00395 


0.00005 


0.4 x 0.24 


2 


13 


0.0048 


0.00002 


0.1 x 0.1 


3 


13 


0.006 


0.001 


4x4 


4 


7 


0.014 


0.002 


0.4 x 0.4 


4 


1.36^ 


0.425 


0.085 


1.8 x 1.3 


5 


0.46 


5.3 


1.8 


2.4 x 1.1 


5 


0.0195 


109.6 


21.9 (b) 


2.3 x 2.3 


6 


0.0117 


5.86 


1.17® 


2.3 x 2.3 


6 


0.0080^ 


3.9 


0.3^ 


2.4 x 2.4 


7 


0.0058^ 


2.63 


0.1^ 


2.4 x 2.4 


7 


0.0047 (6) 


1.8 


0.36 (c) 


2.3 x 2.3 


6 



References. - (1) Marti et al. (1993); (2) Gomez et al. (1995); 
(3) Marti et al. (1999); (4) Gomez et al. (2003); (5) this work; (6) 
Aspin et al. (1994); (7) Qiu et al. (2008). 

(a) The data a t this wavelength were obtained by combining our 
high and low angular resolution data and applying the weighting 
factor, robust=0.5. This image yields a beam area similar to that for 
the image at 456 /zm. 

(6) Points not included in the fit of the SED of MM1. 

( c )The errors were estimated to be 20%. 
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Table 5. MM1 SED fit 



Source 




T456/im P 


T 

A dust 




Luminosity 




n 2 




K 


Jy 




MM1 


1.1±0.4 


0.22±0.08 0.53±0.05 


109±3 


0.0064 ±0.0004 0.18±0.01 


3300±500 



( a ) Free-free emission contribution at 456 /iin in Jy. 



Table 6. Estimated masses and column densities. 



Source 


Beam 


Frequency 




Tdust 


Diameter 


M 




n (b) 






(GHz) 


(cm 2 g 1 ) 


K 


// 


(M ) 


(xlO 23 cm- 2 ) 


(XlO 7 cm- 3 ) 


MM1 


CK'68 x Of/38; 13?0 


220 


0.97 


109 


0.36 


4.1±0.2 


113 


185 


MM1 


8'.'1 x 3'.'0; 34? 1 


220 


0.97 


109 


2.1 


5.6±0.2 


4.7±0.3 


1.33 


MM1 


2" 4 X l'.'l; 37? 1 


658 


1.74 


109 


1.17 


3±1 


9 


5 


MM2(E) 


Of.' 68 x tf'38; 13?0 


220 


0.93 


35 


0.36 


1.9±0.2 


55.0 


90 


MM2(W) 


Of.' 68 x 0'.'38; 13?0 


220 


0.93 


35 


0.7 


1.7±0.3 


13±1 


11 


MM2 


8'.'1 x 3'.'0; 34? 1 


220 


0.93 


35 


2.7 


14.9±0.4 


7.6±0.4 


1.7 


MM2 


2'.'4 X l'.'l; 37? 1 


658 


5.35 


35 


1.17 


17±2 


46.7 


24 



("'Some column densities are lower limits because the unresolved nature of the sources and do not include uncertainties (the same 
applies for number densities). 

( b )lt is assumed a spherical shape of the sources for the number density estimations. 
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Table 7. Properties of the millimeter sources 



Source 




R 


Nh 2 


M 500 (fe) 


M 250 o (c) 


T 


L 








(AU) 


(cm -2 ) 


(M ) 


(M ) 


(K) 


(L ) 


yr 


MM1 


0.5 


< 136 


1.1 x 10 25 


4.1 


5.6 


110 


> 3.3xl0 3 


10 5 


MM2(E) 


1.6 


< 300 


5.5 x 10 24 


1.9 


14.9 


35^ 




2 x 10 3 


MM2(W) 


(/) 


600 


1.3 x 10 24 


1.7 


(/) 


(/) 







( a ) Dust opacity index. 

( b ) Mass of the dust inside a radius of 500 AU (from the high angular resolution 
observations). 

(^Mass of the dust inside a radius of 2500 AU (from the low angular resolution 
observations). 

^The dynamical time of the outflows associated with MM1 and MM2(E) (Qiu & 
Zhang 2009). 

(^Temperature given in Qiu & Zhang 2009. 

^These quantities are estimated for MM2(E) and MM2(W) together. 



